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ON PARAMETER IDENTIFICATION FOR LARGE SPACE STRUCTURES 


By 

S.M. Joshi.1 and G.L. Goglia 2 


ABSTRACT 


The design of a controller for large space structures (LSS) based on 
the LQG theory requires the knowledge of the LSS parameters. Since apriori 
knowledge of the parameters is usually not reliable, the parameters must be 
identified prior to the controller synthesis, using methods such as the 
maximum likelihood technique. In this report, an expression is obtained for 
the Fisher information matrix for LSS, from which Crame'r-Rao bounds can be 
obtained in order to determine the accuracy with which the parameters can be 
identified. 


INTRODUCTION 


The problems encountered in controlling large flexible space structures 
arise mainly because of (1) infinite number of structural modes, (2) numer- 
ous low-frequency structural modes, (3) extremely small structural damping 
ratios, and (4) closely-spaced structural mode frequencies. The controller 
design for large space structures (LSS) can be accomplished by using tools 
of modern control theory. Two types of control systems synthesis techniques 
were discussed in reference 1. The first technique, the "collocated con- 
troller technique," utilizes the stability and robustness properties of 
controllers employing positive definite feedback and collocated actuators 
and sensors. The second technique results in a reduced order Linear-Quadra- 
tic Gaussian (LQG) controller, where the order reduction is accomplished 
by methods such as those described in reference 2. 
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In the course of the present grant, these two techniques were applied 
(ref. 1) for synthesizing control system for a large space antenna, and RMS 
(root-mean-square) performance, described by pointing, surface and feed 
misalignment errors, was obtained in the presence of sensor and actuator 
noise. The results were also reported in the preceding progress report of 
this grant. One of the outcomes of this study was that the LQG controller 
yielded markedly superior performance. However, the design of an LQG 
controller requires at least a reasonable knowledge of the parameters, i.e. 
natural frequencies, damping ratios and mode shapes (or slopes). It is 
extremely difficult, if not impossible, to accurately determine these 
parameters apriori, on the ground. Therefore, the parameter identification 
will probably have to be performed in orbit, and control system designed 
subsequently. 

In the light of the above discussion, it is important to investigate 
the identifiability of the LSS parameters. In this report, expressions for 
the information matrix for LSS are derived using rate gyros as sensors. 
Crame'r- Rao bounds, which indicate the best theoretical accuracy with which 
parameters can be identified, can be obtained by inverting the Fisher 
information matrix. 


FISHER INFORMATION MATRIX FOR LSS 
Consider a linear constant coefficient system described by 

x = Ax + Bu + v (1) 

c 

y * Cx + Du (2) 

z * y + w (3) 

where x, u, y. and z represent n x 1 state vector, m x 1 input vector, 
l x 1 output vector and 9. x 1 observation vector; v c and w c repre- 
sent process and observation noise (stationary, zero mean, white). A, B, C 

and D are functions of a parameter vector p which has the dimension n x 1 

P 


The discrete-time, maximum likelihood identification problem can then be 
stated as follows: 

Given observations z(t) at t ■ kT, k * 0,1,2,. ..,N, for a known in- 
put u(t), determine the estimate $ of p which maximizes the conditional 
density f(Z^/p) where 

Z N ' ‘V *1 V 

If the apriori density of p is known, then the conditional density 
(likelihood) above can be replaced by the unconditional density. 

The discretized version of equations (1), (2), (3) is given by: 


x(k + 

1) - F x(k) + Gu(k) + v(k) 

(4) 

y(k) 

* C x(k) + Du(k) 

(5) 

z(k) 

■ y(k) ♦ w(k) 

(6) 


If x(o), v(k) and w(k) are Gaussian, the density f(Z^/p) is Gaussian. (If 
f(p) is Gaussian, then the unconditional density f(Z^) is also Gaussian.) 

N ! 

f(Z N /p) * ci exp [-1/2 j£q ( z(k) - y(k)) W (z(k) - y(k)l] (7) 

where ci is a constant and W is the covariance matrix of w(k). 

It can be proved that (ref. 3) 

* T -1 

E [ (p-$) (p-p) /p] > J (8) 

where $ is any absolutely unbiased estimate of p, and where 

J ’ e {[3V" £( Vp>] T |)p } '* 

(9) 
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It can be proved from equations (7) and (9) that OftlG INAL P,' 1 "" 
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J is known as the "Fisher information matrix," and J“ 1 is the Cramer-Rao 
lower bound on the covariance of the parameter estimation error. 


App l ication to LSS .- The unforced modal equations of motion for a 
flexible LSS are given by: 


q ' + 6 . + w?q. * 0 

i i M i i M i 


(i ■ 1, 2, 3 n ) 

hi 


(id 


where n is the number of modes. The output equation is: 
m 


1 "iilVi 


( 12 ) 


where q^, P., end ^ are the modal amplitude, damping ratio, natural 


frequency, and (4x1) mode shape (or mode slope) matrix for the ith struc- 
tural mode, y is the 4x1 output vector. (The rigid-body equations are 
not given because they do not play a role in the identification of the 
structural parameters. It is assumed that the mass and inertias of the LSS 
are known). The parameters to be estimated for the ith mode are: 



• • » 



(13) 


Thus there are (1+2) parameters per mode, or a total of n ( 1+2) para- 

n 

meters to be estimated. 

The evaluation of the Fisher information matrix in equation (10) 
requires the computation of ^ . 

Since the parameters to be identified appear explicitly in the 
continuous model, it is more convenient to first write the continuous 
equations for the sensitivity states (9q^/3p), and then to discretize them. 
The sensitivity state equations are (from Eq.ll): 
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(14) 

(15) 


If rate 


gyros are used for rate measurements , 

_ n 


.Z® ♦. 4. 

i*l i ’i 


(16) 


where is the i x 1 mode-slope matrix for the ith mods. (There would 
be three y vectors corresponding to each axis; only one is given here to 
retain simplicity). 3y/3Pj* 3y/3w^, 3y/3$^ (j"l,2,..,A) can be obtained 

from (14) (15) and (16). The evaluation of 3y/ 3p then requires the 
discretization of the coupled continuous differential equations (11), (14) 
and (15), which represents a sixth-order system (for each mode). Therefore 
the total order of the system is 6n . However, since the equations for each 

IB 

mode are not coupled with the other modes, discretization for each mode can 
be carried out separately. 

The sampling period is next selected, so that the sampling frequency is 
higher than the Nyquist frequency based on the highest mode estimated. The 
estimation time (N in Eq. 10) is determined so that the final time is at 
least equal to the largest 1/(P<») (i.e., the longest decay time). Nonzero 
initial conditions are created by aoplying torques through the torquers 
(e.g. on the mast of the hoop/column antenna). The applied torques could 
consist of sine waves, or pulses, etc. The torques are turned off after the 
modes are sufficiently excited (as indicated by the deflections) , and Fisher 
information matrix J is computed as the unforced motion continues. The 
Cramer-Rao standard deviations for each parameter are then the square roots 
of the appropriate diagonal elements of J~*. 


CONCL.JDINC REMARKS 

In order to investigate the accuracy with which LSS parameters can be 
identified, expressions were derived for obtaining the Cramer-Rao bounds. A 
computer program has been written to compute the bounds, and is presently 
being used to investigate the identifiability of the parameters of the 122m 
hoop/column antenna. The results will be reported in a future report. 
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